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I. INTRODUCTION 



One of the outstanding problems for high-energy physics remains the understanding of the dy- 
namics of flavor. Existing experimental results on kaon, D-meson, 5-meson, and hyperon decays, 
as well as neutral-meson mixing, are all consistent with the loop-induced nature of flavor-changing 
neutral currents (FCNC's) in the standard model (SM) and also with the unitarity of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix with three generations. The continuing study of these processes 
with increased precision will play a crucial role in the search for physics beyond the SM. 

In many types of new physics, the new particles are heavier than their SM counterparts and their 
effects can be described by an effective low-energy theory. A complete set of operators of dimension 
six describing deviations from the SM has been presented in Ref. [1]. A less ambitious program is to 
study only those operators that appear when the new physics effectively modifies the SM couplings 
between gauge bosons and certain fermions [2]. The case of anomalous top-quark couplings has 
been treated before in the literature [3, 4], and it was found that they are most tightly constrained 
by the b ^ s'y decay. Interestingly, this mode does not place severe constraints on anomalous 
charm-quark couplings due to the relative smallness of the charm mass. 

In this paper we focus on new physics affecting primarily the charged weak currents involving the 
charm quark. Including the SM term, the effective Lagrangian in the unitary gauge for a general 
parametrization of anomalous interactions of the W boson with an up-type quark [/^ and a down- 
type quark Di can be written as 

^UDw = -^V,,U,r[{l + ^li)PL + 4^r]AW^; + H.c. , (1) 

where g is the weak coupling constant, we have normalized the anomalous couplings Kf^l relative to 
the usual CKM-matrix elements V^;, and -Pl r = K-'-'F^s)- general, Ki are complex and, as such, 
provide new sources of CP violation. In Appendix A we discuss the general parametrization of the 
quark- mixing matrix underlying Eq. (1), paying particular attention to the number of independent 
parameters that are allowed. 

In addition to affecting weak decays though tree- level interactions, the new couplings in Eq. (1) 
modify effective flavor-changing and -conserving couplings at one loop. In this work, we evaluate 
several one-loop transitions induced by the new couplings via magnetic-dipole, penguin, and box 
diagrams. We include the operators generated this way in our phenomenological analysis. 

We present a comprehensive picture of existing constraints which shows that deviations from 
SM couplings at the percent level are still possible, particularly for right-handed interactions. Our 
study will also serve as a guide as to which future measurements provide the most sensitive tests for 
new physics that can be parametrized with anomalous IV-boson couplings to the charm quark. For 
the CP- violating (imaginary) parts of the couplings, the electric dipole moment of the neutron and 
the hyperon asymmetry A^^^ are the most promising channels to probe for right-handed couplings, 
whereas more precise measurements of sin(2/3) and sin(2/3s) are the most promising probes for 
left-handed couplings. Constraints on the real parts of the right-handed couplings can be further 
improved with better measurements of semileptonic B and D decays. 
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II. ONE-LOOP PROCESSES 



In this section, we collect the main formulas for the loop-induced processes we are considering. 
All our calculations are performed in the unitary gauge, and the results have been compared to 
existing ones, where available. We have summarized our loop calculations in Appendix B. 

A. Dipole penguin operators 

Of particular interest are the electromagnetic and chromomagnetic dipole operators, which can 
give rise to potentially large corrections to SM processes [5, 6] and be expressed as 

= Y^2 id\uPKd ± d\^PLd) , (2) 

Qt = ^2 id\utaPKd ± d\^t^P^d) Gf , (3) 

where d and d' ^ d are down-type quarks, F^^ and are the usual photon and gluon field- 
strength tensors, respectively, with e and being their coupling constants, and Tr(tatfe) = \5ab- 
The main effect of these operators is generated at one loop with the W boson coupling to the left- 
handed current at one vertex and to the right-handed current at the other vertex. It corresponds 
to the terms linear in the k^'s in Eqs. (Bl) and (B2), leading to the effective Hamiltonian 

H,,, = c+g+ + C7g; + c;g+ + c;g; + h.c. , (4) 

where at the H^-mass scale 

q=u,c,t 
q=u,c,t 

with Gp being the Fermi constant, = 'rhg{mg^ /^w^ ^^'^ given in Eqs. (B6) and (B7). 

The contributions of these operators are potentially enlarged relative to the corresponding ones in 
the SM due to the enhancement factors of m^/m^ and m^/rrif^ in s ^ d and b {d, s) transitions, 
respectively, and also due to Fq^x^) and GQ(Xg) being larger than their SM counterparts. For the 
case of anomalous tbW couplings, the formulas above agree with those found in the literature [3] . 

B. Electric dipole moments 

The flavor-conserving counterparts of Q^^g above contribute to the electric and color dipole- 
moments of the d and s quarks. Based on Eqs. (4) and (5), one can write the effective Hamiltonian 
for such contributions to the dipole moments of the d quark as 

nr = hdf-da^^^.dF^'^ + ^V^aTs^^r , (6) 
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where in our case 



^ q=u,c,t 



q=u,c,t 



at the m^r scale. These expressions agree with those derived from quark-l^ loop diagrams in a left- 
right model [7], after appropriate changes are made. The corresponding quantities for the s quark 
are similar in form. Since the dipole moments ^ contribute to the electric dipole moment of the 
neutron [8-10], it can be used to place constraints on ImK^^^. 



C. Operators generated by Z-penguin, 7-penguin, and box diagrams 



The anomalous quark-W^ couplings also generate flavor-changing neutral-current interactions via 
Z-penguin, 7-penguin, and box diagrams. They will therefore affect other loop-generated processes, 
such a.s K —>■ TTvv, Kl — , and neutral-meson mixing. 

The effective theory with anomalous couplings is not renormalizable, and this results in diver- 
gent contributions to some of the processes we consider. These divergences are understood in the 
context of effective field theories as contributions to the coefficients of higher-dimension operators. 
These operators then enter the calculation as additional 'anomalous couplings', introducing new 
parameters to be extracted from experiment. For our numerical analysis, we will limit ourselves to 
the anomalous couplings of Eq. (1), ignoring the higher-dimension operators. In so doing, we trade 
the possibility of obtaining precise predictions in specific models for order-of-magnitude estimates 
of the effects of new physics parametrized in a model-independent way. We will rely on the common 
procedure [11] of using dimensional regularization, dropping the resulting pole in four dimensions, 
and identifying the renormalization scale /x with the scale of the new physics underlying the effective 
theory. Our results will thus contain a logarithmic term of the form / m^r) in which we set 
/i = A = 1 TeV for definiteness. In addition to the logarithmic term representing the new-physics 
contribution, we have also kept in our estimates those finite terms that correspond to contributions 
from SM quarks in the loops. ^ 

We consider the contributions of the anomalous couplings to dd' — * uu, dd' and 
dd' dd', relegating the main results of the calculation to Appendix B. It follows that the 
effective Hamiltonians generated by the anomalous charm couplings are at the m^r scale 



a 



+ 



TT sin^ 



(-3 \n-^ + 4X,{x,))dYP^du^^P, 
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v/Svr sin^e. 
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(4x^ - 3) In 
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W 



dYPndiyi.PL'y 



(8) 



^ Explicit examples of the type of divergence cancelation resulting in a logarithmic term as described above can be 
found in Ref. [12], where the new-physics scale is given by the masses of non-SM Higgs-bosons. 
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(10) 

where li' 7^ li, we have kept terms hnear in and quadratic in k^, A^ = Vg^/V^^, and is the 
Weinberg angle. The functions Xq, Yq, Zq, X, Y, Z, and ;Bi 2,3 can be found in Appendix B. 



III. TREE-LEVEL CONSTRAINTS 



From now on, we focus on the anomalous charm couplings k^^ neglecting the corresponding 
u and t anomalous couplings. To obtain constraints on the couplings, we begin by exploring their 
tree-level contributions to three different sets of processes, -D(s) iu, exclusive and inclusive 
b — >■ ci'u transitions, and mixing-induced CP violation in B — > J/ipK and B Vc-^^^ where the 
couplings may play some interesting roles. 



A. D, Ds iv 



From the Lagrangian in Eq. (1), at tree level one derives the decay rate 

TiD^iu) = 9l3^li^(i_l4 \ \VUl + -^cd--fd)\\ (11) 

where the decay constant is defined by {0\d'-f^^c\D{p)) = ifj^p^. Changing V^^^ i\,d^^ 
fjj to V^g, ^cs^, rriD^ and /d^, respectively, one obtains the decay width T{Dg — > iu). 
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Recent measurements of Dg iv yield [13, 14] 

/^"P = (205.8 ± 8.9) MeV , /^"^ = (261.2 ± 6.9) MeV , (12) 

whereas SM calculations give [15, 16] 

f% = (202 ±8) MeV, fj^^ = (240 ± 7) MeV . (13) 

Evidently, for D ^ iu the data agree with theoretical predictions well, but for Dg iu there is 
deviation at the 2-sigma level. It has been argued that this deviation may be due to physics beyond 
the SM [15], but it is too early to conclude that new physics is needed. 

Nevertheless, one can turn the argument around to constrain new physics by assuming that 
the discrepancy between the calculated and measured values of the decay constants arose from 
the anomalous couplings, as the r(D(s) — > iu) formulas would imply. Using the experimental and 
theoretical numbers above, one can then extract 

\Re{K];,-^f,)\ < 0.04, (14) 
< Re(«:L < q.I . (15) 

B. Semileptonic B decay and extraction of Vet 

The interaction in Eq. (1) will also affect the extraction of Vcb from semileptonic B decay. At the 
quark level, the effect of the new couplings is to scale the hadronic vector and axial-vector currents 
by the factors 1 + ± k^, respectively. This has the following implications. 

First, the semileptonic exclusive decay B Dev^ is sensitive only to the vector form-factor, 
and thus the differential (and total) decay rate simply gets multiplied by 1 1 + kJiJ, + k^]^. To linear 
order in the k's, this means that what is measured in this mode is 

= (1 + Re k]:, + Re kI) = (39.4 ± 4.4) x lO'^* . (16) 

The number above and the other ones below for are quoted from Ref. [17], and their errors 
result from adding the experimental and theoretical uncertainties given therein in quadrature. 

Second, the semileptonic exclusive decay B D*eu^ is sensitive to both the vector and axial- 
vector currents. In the heavy-quark limit, w = v ■ v' = 1 (where v and v' are the four- velocities of 
the B and D*, respectively), only the axial- vector current survives [18], and so the decay rate in this 
limit would simply get multiplied by |l + — One can do better than this by considering 

the form factors in the heavy-quark effective theory (HQET) where they either vanish or can be 
written in terms of the Isgur-Wise function C.i'w) with the normalization ^(1) = 1 [18]. Treating 
the form factors as constants throughout the kinematically allowed range 1 < w < 1.5, one then 
finds to linear order in k 

Vf = V;fe(l + Refi:)i-0.93ReK)|) = (38.6 ± 1.4) x 10"^ . (17) 
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Third, the semileptonic inclusive decay rate can be easily calculated to be 

T{b^ce~u^) = ^|K5p{i^W(|l + «:^,r + hf,r)+2G(r)Re[(l + «:^,)«:f;]} , (18) 

where r = m^/m^ ^ 0.3, 

F(r) = 1 - Sr^ + 8r^ - - 24r^ Inr , 

G{r) = -8r [l + 9r^ - 9r^ - + 12r^ (l + r^) In r] . (19) 
It follows that to linear order in k, 

= (l + Re/tJi - l.bRen^) = (41.6 ± 0.6) x 10^^ . (20) 

From these results it is evident that it is not possible to extract a bound on kJj^ (as long as 
quadratic effects are ignored), but we can extract bounds on k^. For example, we can do a two- 
parameter fit to Eqs. (16), (17), and (20) to find a minimum for 

V^f, {1 + Re K^^) = 0.038, ReK^^ = -0.057, (21) 

with a corresponding 68%-C.L. interval (1-cr error) 

-0.13 < Renf^ < . (22) 



C. CP violation in B ^ J/ipK and B r}^K 

One of the decay modes expected to provide a clean determination of the unitarity-triangle 
parameter (3 from the measurement of time-dependent CP violation is -B Vc-^^ just like B — > 
J/ipK. The SM predicts the same sin(2/3) for the two processes, whereas the current data for its 
effective values are [19] 

sin(2/?;^) = 0.657 ±0.025, sin(2/3^y = 0.93 ± 0.17 , (23) 

which disagree with each other at the 1.5-sigma level. Once again we can use the difference between 
the two measurements to constrain the new physics parametrized by the anomalous couplings. Since 
sin (2/3*^^) measures the difference between the phase of the S-mixing matrix element and the phase 
of the ratio of amplitudes for the B decay and its antiparticle decay [20], then the discrepancy in 
j3'^^ between the two modes must arise from a difference between the phases of their amplitude 
ratios. 

The effective Hamiltonian for the b —>■ sec transition including the contribution of anomalous 
couplings can be written as 

+ Cf L -Sr^YP^C^ Cnl.PlPm + Cf"^ sYPkC .Plp) , (24) 
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where (j^^'^^") g^^e the Wilson coefficients, m and n are color indices, and we have neglected penguin 
operators. To linear order in k, the Wilson coefficients at the scale are C2{^m^r) = 1 + n^* + n^f,, 
C2^(m^) = K^, and Cf■^{m^r) = nf*. These can be evolved down to a renormalization scale 
fi ~ to become [21, 22] 



(25) 



at leading order in QCD, where r] = a^{m^^)/a^{fi), 

To determine the amplitudes for B ^ J / ipK, rj^K, we adopt the naive factorization approxima- 
tion. The relevant matrix elements and parameter values are collected in Appendix C. It follows 
that 



where 



(26) 



^ ^LR,RL 
^ _ n I ^LR,RL _ ^LR,RL , '^2 /07^ 

The presence of the second term in the B — > rj^K amplitude offers the possibility of sm(2P'^^^ in 
this decay mode being different from that in i? — J/ipK. Defining 

(a\^(rL) - af^(rL))ml 

we then obtain to ffist order in k 

Pf^K = /?;1, + arg(l + rJ ^ pf^ + lmr,. (29) 

Taking /i = = 4.2 GeV and A^, = 3, we find a^{n) = 0.076 and r^(/i) ^ 20(^^6 - k^/). Since 
the experimental numbers in Eq. (23) imply 

pf^^ = 0.60 ± 0.23 , = 0.358 ± 0.017 , (30) 

in view of Eq. (29) we can then impose —0.005 < Im r^(/i) < 0.4, which leads to the bound 
-2.5 X 10-^ < Im(K^ft + K^^) < 0.02. 

It is well known, however, that this naive factorization procedure fails to reproduce the exper- 
imental branching ratios, which can be better fit with N^. — 2 [23]. Using this value we obtain 
instead r^(/i) ~ 8(k^ — k^*). This would increase the upper bound for Im(K^ -|- k^) above by 
about a factor of two, within the intrinsic uncertainty of our calculation, 

-5x10"^ < Im(K|| + K^J < 0.04. (31) 
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IV. CONSTRAINTS FROM DIPOLE PENGUIN OPERATORS 



We turn next to constraints from the magnetic-penguin transitions d — > (i'7 and d — d'g. The 
specific processes we discuss are b — ^ 57, d — > S7, the CP-violation parameters e and e' in kaon 
mixing and decay, and hyperon CP violation. 



A. b 



Including the SM contribution, the effective Hamiltonian for h ^ is 



n 



jjLU 1 



where to 0{n) 



n 



n 



' qs ' qb 



v*v. 

' qs ' qb 



q=u,c,t 



(1 + + m, FqS^^ (x J + /^J Fo (x J 



(32) 



(33) 



following from Eqs. (Bl) and (B3). The corresponding expressions for b ^ sg are similar in form 
and follow from Eqs. (B2) and (B4). 

The experimental data on b s'j have been found to impose very strong constraints on k^^^^, 
limiting them to below the few-percent level [3]. Since V^^Kb — ~^t*s^tb ''^t ^ 
preceding equations indicate that, if all the k's were comparable in size, the top contributions 
would be larger than the charm ones by almost two orders of magnitude. All this means that 
6 ^ 57 offers relatively weak constraints on k^^^^, with upper bounds at the level of 0{1). 



B. s ^ d-y 



In an analogous manner, the anomalous couplings contribute to the short-distance transition 
s — d'-f, but in this case the charm contribution is expected to be more important than the top 
one. At lower energies, and Cg mix because of QCD corrections. At yU = 1 GeV we have [22] 



where 



V 



2/23 / \2/25 / \2/27 



Numerically, keeping only the charm contributions yields 



(34) 



(35) 



-38 + 0.023 i)(K^,±K^;) X lO^^GeV"^ . 



(36) 



Hyperon and kaon radiative-weak decays provide the relevant constraints, the former being 
somewhat stronger and yielding [24] 



87r2 Gp 



< 12MeV . (37) 



This translates into 

I + I ^ 3 , (3^ 
which is a very weak bound compared to Eqs. (14) and (15). 



C. e and e' 



The gluonic dipole operators contribute to the CP- violation parameters e and e' in kaon mix- 
ing and decay, respectively. Since is parity conserving, it contributes to e via long-distance 
effects [25, 26]. Being parity violating, Q~ contributes to e'. One finds [6, 26, 27] 

(e)^ = -2.3 X 10^GeVSJmC+(/i) , (^^^ = 4.4 x 10^ GeV B,, Im C-(/i) , (39) 

where the contributions of the anomalous charm couplings to are 

C^i^,) = rfCfimy,) = (-21 + 0.013 z) («:^, ± K^;) X lO^^GeV^i (40) 

for /i = 1 GeV and the hadronic uncertainties are represented by 

0.2 < |BJ < 1 , 0.5 < \B^,\ < 2 . (41) 

The experimental data are |e| = (2.229 ± 0.012) x 10"^ and Re(e7e) = (1.65 ± 0.26) x 10"^ [13]. 
The SM predicts |e|sM = {^mt^^^^) x 10"^ [28], but for e' the SM calculation still involves a large 
uncertainty [29]. Consequently, we require that 



|e|^ < 0.7xlO"S j^^J < 1.7x10"^. (42) 

The resulting constraints on k^^^ are complicated and will be presented in Fig. 2 in Sec. VII. 

There are other loop-generated operators contributing to e, and hence they provide more con- 
straints on the anomalous charm couplings. These operators will be discussed in Sec. VI C 



D. Hyperon nonleptonic decays 

Hyperon decays provide an additional environment to study CP-violating | AS"! = 1 interactions. 
The main observable of interest in this case is the CP- violating asymmetry A= (a + a) / (a — a) , 
where a is a decay parameter in the decay of a hyperon into another baryon and a spinless meson 
and a is the corresponding parameter in the antiparticle process [30]. 
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Experimentally, a preliminary value A^^i = (—6 ± 3) x 10 ^ for A measured in the decay chain 
S — > Air — pirir has recently been reported [31]. The SM prediction is < 5 x 10"^ [30], 

which is an order of magnitude smaller than the central value of the measurement. Since this is 
only a 2-sigma disagreement, it is premature to attribute it to new physics. However, this difference 
can also be used to constrain the anomalous charm couplings. 

The contribution of the gluonic dipole operators to the asymmetry A^^ has been estimated in 
Ref. [27]. The result can be written as 

(^-a). = 10^5+lmC;(/i) + 10^i?„lmC;(/i) , (43) 
where C^(/i) due to the anomalous charm couplings are given in Eq. (40) and 

-1.4 < < 0.5 , -0.9 < 5„ < 1.3 

reflect the hadronic uncertainties. The preliminary HyperCP result above suggests that 

-9 X 10"^ < (A=a)^ < -3 X 10"^ . (45) 
The resulting constraints are shown in Fig. 2 in Sec. VII. 



(44) 



V. ELECTRIC DIPOLE MOMENT OF NEUTRON 



The flavor- conserving counterparts of the magnetic-dipole operators discussed above contribute 
to the neutron EDM. The latter is described by the effective Lagrangian 



(46) 



The dipole moments d^'^™ and d^'^™ of the d quark in Eq. (7) contribute to d^. Using the valence 
quark model, we have [8] 



4'^ = ldf-{li) + tedf^{ii) , 

where at /i = 1 GeV 

^^(/i) = r/«rff-(m^) + l{rf -rf)edf^{m^) , 
df^{^^) = fdf^im^) , 
with f] being given in Eq. (35). The anomalous charm contribution is then 

d^f) = 6.9ImK|;^^ X 10^22 gpj^ _ 



(47) 



(4J 



(49) 



Similarly, the electric and color dipole-moments of the s quark produced by the anomalous charm 
couplings are 



rff-(/i) 



-eGr: 



K„ Ini K^, 



(50) 



2^2 7r2' 
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and so their contribution to is given by 

rfW = B^df''\fi) + B^edf"'{fi) = (82 5, + 46 5jlmfi:^, X 10-22ecm , (51) 

(s) 

where —0.35 < < —0.01 and 0.01 < B^ < 0.26 reflect the wide range of estimates for dn in 
the hterature [10], in contrast to those for dlfK The combined contribution of dlt'^^ is then 

K), = 4'^ = (0.69Im«:^, + BJm^l^ x lO^^igcm , (52) 

where 

-2.8 < 5„ < +1.1 . (53) 
Prom the experimental bound |(i„|exp < 2.9 x 10~^^ ecm at 90% C.L. [13], we will impose the bound 

\dj^ < 2.9 X 10-2^ ecm (54) 

in Sec. Vll to restrict the anomalous couplings further. 

The s-quark dipole moments d^'^^ and c/^*^™ above also contribute to the EDM of the A hyperon 
and therefore may be constrained directly by experiment. However, the experimental limit, d^^ = 
(-3.0 ±7.4) X 10"^^ ecm [32] or d^ < 1.5 x lO'^^ecm at 90% C.L. [13], is very weak compared to 
Eq. (54) and hence will not be used for constraining the couplings. 

VI. OTHER LOOP CONSTRAINTS 

In this section we explore several other processes where the anomalous charm couplings can 
contribute via penguin and box diagrams. 



A. K+ 7V+ 



TV ' UU 



To quantify the contribution of the anomalous charm couplings to this mode, it is convenient to 
compare it with the dominant contribution in the SM. The latter comes from the top loop and is 
given by [21, 33] 

-^sm(A'^-^M = sin^ 9^ ^"^"^^^ {'n^\si ,d\K^) vi^{\ - 75)^ , (55) 

following from Eq. (B9), without QCD corrections. It is also convenient to neglect the masses of 
the leptons associated with the neutrinos in the new contribution, as in Eq. (8), so that we can 
work with just one of them. The total amplitude can thus be written in terms of the SM amplitude 

above as 

M{K^^T^^vv) = (1 + 5) A^sm(^^ ^ TT+Z/I^) , (56) 
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where to linear order in k 

c ^ VcdV^s (^cd + ^Ys) [-3 ln(A/m^) + AX,{x,)] 

In the above expression we have used ^(xj ~ 1.4 instead of Xq(xj) in the denominator to 
incorporate the QCD corrections in the SM [34]. The SM prediction for the branching ratio is 
B's^y^iK^ n^uu) = (8.5 ± 0.7) x 10"^^ [34], to be compared with its experimental value B^^^ = 
(l-73ll:o5) X 10"^° [35]. Accordingly, we require —0.2 < Re5 < 1, which translates into 

-2.5xl0~^ < -Re{K^a + k]:^) + 0A2lm{K^^- k]:^) < 1.3 x lO'^ . (58) 

It is interesting to compare the anomalous charm contribution to the SM charm contribution. 
Their ratio is 

■M,(j^+-.7r+z/P) _ {k^^, + K^:) [-3 ln(A/m^) + AX,{x,)] ^^^^ 
M^il{K+ ^ Ti+uu) 

where 6 x 10~^ ^ ~ 1 ^ 10~^ incorporates QCD corrections and lepton-mass dependence [21]. 
With the k];^ + hi]:^ at the upper end of the range above, the two contributions are similar in 
size. This implies that the current experimental situation admits a 100% uncertainty in the charm 
contribution to the branching ratio, much larger than the theoretical uncertainty within the SM. 

Constraints on the anomalous couplings can also be extracted from the related 5-meson modes, 
B Xuu, but the resulting bounds are about three orders of magnitude weaker due to unfavorable 
CKM angles. Experimentally, only upper limits for their decay rates are currently available [19]. 
For these reasons, we do not discuss them further. 



4X, 



NL 



B. Kl 

The dominant part of the short- distance contribution to the SM amplitude for — > fi'^fi~ is 
again induced by the top loop and can be expressed as [21, 33] 

Ml^{K' ^ /iVi = 2^ V,,VZ Y,{x,) {Q\-sri,d\K') , (60) 

from Eq. (9). Combining this with the anomalous charm contribution in Eq. (9), we arrive at the 
total short-distance amplitude 

-Msd(^l^/^V) = (l + 5')-MiM(^L^/^V) , (61) 

where to linear order in /t 

^, ^ Re[K*As (^L + ^cd*)] [-3 ln(A/m^) + 4^0 (x,)] 

ARe{V:,V,,)Y{x,) 

with Y[x^ ~ 0.95 being the QCD-corrected value of ^o(^t) [36]- Since the measured branching 
ratio, B{^K^ — > fj.^fj.~) = (6.84 ±0.11) x 10"^ [13], is almost saturated by the absorptive part of the 
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long-distance contribution, B^-^^^ = (6.64 ± 0.07) x 10~^ [37], the difference between them suggests 
the allowed room for new physics, Bj^p ^ 3.8 x 10"^°, the upper bound being about one half of 
the SM short- distance contribution, B^^ = (7.9 ±1.2) x 10~^° [36]. Consequently, we demand 
\6'\ < 0.2, which implies 

+ + 6 X 10-^ lm{Kl - /^^ | < 1.5 x 10"^ . (63) 

One could also carry out a similar analysis as above for B , but the CKM angles in that 

case are such that the constraints would be much weaker. In addition, only experimental bounds 
on the rates are currently available [19]. 



C. K-K mixing 

The matrix element Mj^g K^-K^ mixing is defined by [21] 

2m ^Mg = {K'\n,,_^,,\K') (64) 

where the effective Hamiltonian Ti^g^^g consists of SM and new-physics terms. For the latter, the 
contribution of the anomalous charm couplings can be derived from Eq. (10), where the last line is 
negligible compared to the second because of the smallness of x^. Thus 



^K,K _ '^F^W f2 yds 



= -(0.090 + 0.03H) ps'^ + + (2.1 + 0.58i)ps-^K^;/t^, , (65) 

where A^^ = V*^Vgg, we have included QCD-correction factors at leading order with f/ given in 
Eq. (35), m^{fi) + mg{fi) = 142 MeV at fi = 1 GeV, and the other parameters can be found 
in Appendix C. Evidently, the inclusion of the second term in M^''^, albeit quadratic in k^, is 
important as it receives large chiral and QCD enhancement with respect to the first term. 

Now, the difference AMj^ between the Kj^ and Kg masses is related to Mj^ = M^/^ + M^^^ 
by AM^ = 2ReM]^ + AM^^, the long-distance term AM^° being sizable [21]. Since the LD part 
suffers from significant uncertainties, we constrain the anomalous couplings by requiring that their 
contribution to AM^ be less than the largest SM contribution, which comes from the charm loop 
and is given by 

^K,SM ^ ^^^^^^ , ^^^^ ^gg^ 

with the parameter values in Appendix C. The result is 
|0.043Re(K^, + Ky +0.015 Im(«:|r,-Ky - Re(4*/^a + 0.28 Im(«:^;«:^J | < 8.5 x 10"^ (67) 
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A complementary constraint on the couplings can be obtained from the CP-violation parameter e. 
Its magnitude is related to by [21] 



|lmM-^|| 



(68) 



where AM'j^^ = (3.483 ± 0.006) x 10"^^ GeV [13] and the small term containing the CP-violating 
phase in the K nn amplitude has been dropped. Since |e|cxp = (2.229 ± 0.012) x 10"'^ [13] and 
|e|sM = (2-06;!:q;53) ^ [28]; we again demand |e|^ < 0.7 x 10"^ for the contribution in Eq. (65). 
This translates into 

|0.015Re(«:^^ + ^y +0.043 Im(«:;r,-«:a - 0.28Re{^f: k^) - lm{f,f: < 2.5 x 10-V (69) 



D. Bd-Bd mixing 

In the SM the matrix element Mj^g ^d~-^d mixing is dominated by the top-loop contribution 
and given by [21] 

/ ROl-iySM _ I rO\ r'2 ^2 

<r = ^ ^' - ^ flm,^ ( V.V.;)= S, (x,) (70) 

following from Eq. (B18), where a QCD-correction factor rj^ has been included and the parameter 
values can be found in Appendix C. In contrast, since V*^V^f^ and V^^^b are comparable in size, the 
anomalous couplings of charm and top may produce similar effects on Mf2, as can be inferred from 
Eq. (10). The anomalous top contributions to B mixing having been studied before [38], we switch 
them off and get, to linear order in k, 

= ^/B.-B.^.5B.AfK, + .a(-Afx.ln^-$^ 
= (1.6 + 0.63 2)ps-i(«:|r, + «:a ' (71) 



where X"^ = V^^K^b have neglected the terms quadratic in because their quark operators 

do not get as much chiral and QCD enhancement as those in the kaon-mixing case. 

The difference AM^ between the masses of the heavy and light mass-eigenstates is related 
to = Mff^ + Mfi"" by AM^ = 2|Mf2| [21]. The measured value AMf^ = (0.507 ± 

0.005) ps-i [13] agrees with the SM prediction, AM™ = (0.563t[!:[]?^) ps'^ [28]. In the presence 
of the anomalous couplings, these numbers are related by 

AMr = AM^ |1 + S,\ , 6, = . (72) 

Accordingly, we impose —0.2 <Re6^< +0.02, which leads to 

-0.031 < Re(KL+/s:y + OAlm{K];, - k]:,) < 0.003. (73) 
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An additional constraint can be extracted from the (5 measurement m B ^ J / ipK. The anoma- 
lous couplings enter (3'^^ via both the mixing and decay amplitudes. Since the mixing parameters 
Pb^ and qb^ are related to Mf^ by QBa/PBa — ^12/1^^21 PO], we have 



From the decay amplitude in Eq. (26), we derive 

having incorporated the J'i'- mixing factor qx/PK = KdKs/ (KdKs) [20]- It follows that 

Qb, MjB'^^Ks) 
Pb, M{B^-^^Ks) 



Upon comparing the experimental value 2(3°^ = 2/3^^ = 0.717 ± 0.033 from Eq. (30) to the SM 
prediction 2(3^^ = 0.753;[|;^| [28], we then require -0.01 < 2 lm{K^^^ - ) - Im(5^ < 0.08, which 
implies 

-1.5x10"^ < 0.4Re(fi;Ji + Ky -0.69 ImK^Tb + ImKlTrf -0.31 ImK^r, < 0.012. (77) 



E. Bg-Bs mixing 

The SM part of the matrix element Mfg for B^-B^ mixing is also dominated by the top contri- 
bution [21], 

^^M ^ ^2^^^^ ^^^^^ ^y^^y.^ 2 ^^^^ _ ^^g^ 

For the anomalous couplings, again the charm contribution alone is 

^iT = fk^B. VbBb. Af i^]:, + ( - Af a:, In - 5^ \f B, (x,, 



247r2 

(53 - 0.95z)ps"i(/t^^ + ;,L*^) ^ ^73) 



where Af = V^^V^,. 

Similarly to the B^ case, we have here 



AMr = AMf^|l + <5 I , 5, = . (80) 



Ml, 



The experimental value, AM|^p = (17.77±0.12) ps ^ [13], is in agreement with the SM prediction, 
AMf^ = (I7.6l};^) ps-i [28]. These numbers allow us to require -0.09 <Re6^< 0.1, leading to 

-0.014 < Re{K^, + K^h) + 0.018 Im(fi:J;', - < 0.015. (81) 
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A complementary constraint is provided by the parameter (3^ in decay, analogously to (3 in 
decay. In this case, the mode of interest is 5° J /'4"t'i which proceeds from the same h sec 
transition as B^ J/ipK. For the mixing factor, we have 



9s. 



_ Mini + 61) ^ (82) 

and for the ratio of decay amplitudes, as in Eq. (75), 

where {J/ip(j))f is one of the CP eigenstates of the J/ip(f) final-state and 77^ its CP eigenvalue. It 
follows that 



The SM yields 2/5f^ = 0.03614+°;°°^^^ [28], but the measurements of B^ J/tl^ct) yield the average 

oeff 



value 2(3f = 2pf, = O.rrtoH or 2.36+|]:3? [19]- I* is again too early to attribute this difference to 



new physics, but it can be used to impose the bound -0.003 < 2 Im(fi:Ji - kJtJ - Im5, < 0.4, which 
yields 

-0.09 < 0.026 Re(K^^ + + Im(fi;^^-fi:y < 7x10"^. (85) 



VII. SUMMARY AND CONCLUSIONS 



We have explored the phenomenological consequences of anomalous W^-boson couplings to the 
charm quark in a comprehensive way. Most of the constraints we have obtained are summarized in 
Table I. In writing them, we have followed the discussion in Appendix A about the independent 
parameters in the quark-mixing matrices and chosen argK^;'^ = 0. Consequently, we have used the 
condition Imnl;^ = in all the results. All the constraints in this table are quoted as l-cx errors, 
but in some cases the theoretical error is only an order of magnitude and this is not reflected in the 
quoted range. The discussion in the text makes it clear whenever this happens. We leave out from 
the table the processes 6 — 37 and s — > ^7 since the resulting bounds are not competitive with 
the rest. 

In Fig. 1 we show the parameter space of the real and imaginary parts of k,]^^ and assuming 
that only one of them is nonzero at a time. This figure indicates that currently the phase of is 
only loosely constrained and ranges from —90° to 90°. In contrast, the phase of n^^^ is unconstrained 
if its magnitude is small (at the 10~^ level). However, larger values of \n]^f,\, at the few percent level, 
are also allowed provided its phase lies in a range roughly between —150° and —56°. 

We treat the constraints arising from the contributions of magnetic-dipole operators to CP- 
violating observables separately and display those in Fig. 2. These observables receive contributions 
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Process 




Constraint 


# 


U — > lU 


(14) 




|Rc(«:^,-kR,)|<0.04 


1 
1 


1 ) , 0, , 

Us — 


(15) 




0<Rc(«;L _^R) <0.1 




— > ctv 


(22) 




-0.13 < Rck)| < 




Q 
O 


JD — > J /%Jl\^rj^J\ 


[61) 




-5 X 10-4 < Im(K^f, + k)?,) 


< 0.04 


A 

4 


— > Tr^uu 


(bti) 


-1.3 X 10-3 < Re(KL^ + KL) +0.42ImKLs < 2.5 x 10^4 







[b6) 




Re(K^s + i^^d) + 6 X 10-4lmK^s 


< 1.5 X 10-4 





A ]\'T 


(57) 


0.043 Re(KL^ + ) _ 0.015ImKL _ Re[Kf*Kf^) -t 


- 0.28Im(Kf^*K^J < 8.5 x 10-4 


7 
( 


e (mixing) 




0.015 Re(K5^^ + K^^) + 0.043 ImK^s - 0.28Re(K^^*K 


I) -lm{Kf*Kf,)\< 2.5x10-^ 


Q 
O 


AMd 


(73) 




-0.031 < Rc{k^^ + kJtJ + 0.4ImKL < 0.003 


9 


sin(2/3) (mixing) 


(77) 


-1.5 X 10" 


< 0.4Rc(K^b + K^^) - 0.69 ImK^^ - O.SllniK^, < 0.012 


10 


AM, 


(81) 


-0.014 < Rc(kL + ) + 0.018 Im(KL _ ) < 0.015 


11 


sin (2 /3s) (mixing) 


(85) 


-0.09 < 0.026 Rc(kL + ^ ^ lni{K^^ 


- ) < 7 X 10-4 


12 



TABLE I: Summary of constraints, with their equation numbers, from various processes. 



from the anomalous couplings Im/t^^^ that are much larger than SM contributions to the dipole 
operators. These enhanced contributions to e, e', the neutron EDM, and have been studied 
before as they arise within LR models and supersymmetry [6, 26, 39]. 

The calculations for all of these CP-violating observables suffer from large theoretical uncer- 
tainties which we have parametrized with B factors in this paper. For illustration, we display two 
plots in Fig. 2 resulting from choosing two representative sample sets of values of the parameters 
-Be_e'_+^_,n withiu their ranges in Eqs. (41), (44), and (53) and imposing the constraints in Eqs. (42), 
(45), and (54). In each of the plots, the very lightly shaded (yellow) band satisfies the e'/e con- 
straint, the lightly shaded (pink) band the e constraint, the medium shaded (green) band the 
constraint, the heavily shaded (blue) band the constraint, and the dark (red) region all of the 
constraints. It is worth noting that there is a significant amount of the parameter space where 
all of the constraints can be simultaneously satisfied and that the values of Im involved are 
typically of order a few times 10~^ or less. Furthermore, as is obvious from the plots, the neutron- 
EDM constraint is the most restrictive. Also, interestingly, the allowed region of parameter space 
easily accommodates an much larger than the SM prediction, as hinted at by the preliminary 
measurement by HyperCP [31]. 

In order to gain some insight into the constraints in Table I and Fig. 2, we have extracted the 
ranges corresponding to taking only one anomalous coupling at a time to be nonzero (and only for 
the cases of a purely real or a purely imaginary coupling). They are collected in Table II. This 
table shows that, in general, the left-handed couplings are much more constrained than the right- 
handed couplings. Similarly, the imaginary part of the couplings is more tightly constrained than 
the corresponding real part. The largest deviations allowed by current data appear in the real part 



18 




-0JCB02 -0.«B] OjOOW OjCBO] OjO(B2 -0j04 -OjM OjOO OjCH OjM 

FIG. 1: Parameter space of the real and imaginary parts of k^^^ and subject to the relevant constraints 
in Table I, under the assumption that only one k is nonzero at at time. The heavily (blue), medium 
(green), and lightly (yellow) shaded areas in the left plot satisfy constraints #2, #6, and #8, respectively. 
The heavily (blue), medium (green), and lightly (yellow) shaded in the right plot satisfy constraints #9, 
#10, and #12, respectively. The dark (red) region in each plot satisfies all the constraints in it. 



< ReK^^ < 1.5 X 10-^ 


#6 






< ReK^, < 1.5 X 10""^ 


#2, #6 


-6 X 10"^ < Im k);', < 6 X 10"^ 


#8 


-4 X 10~3 <ReK^^<3x IQ-^ 


#9, #10 


-0.02 <lmK^f^<7 X lO"'' 


#10, #12 


-0.04 < ReK^^ < 0.04 


#1 


-2 X IQ-^ < Imnf^ < 2 x 10"^ 


Fig. 2 


-0.1 < Renf^ < 


#2 


-5 X 10""^ < ImK^, < 2 X 10~3 


#4, Fig. 2 


-0.13 < RckJI < 


#3 


-5 X 10"^ < ImK^j < 0.04 


#4 



TABLE II: Constraints on each of the anomalous charm couplings, extracted from Table I and Fig. 2. 

of the right-handed couplings, which can be as large as 10% of the corresponding SM couplings. 

For specific model building, it is useful to recall that the new physics in the quark-mixing matrices 
is being parametrized here as the product of the anomalous coupling and the corresponding CKM 
matrix element, as can be seen in Eqs. (A2) and (A3). The allowed new physics in the left-handed 
sector is at most of (9(A^) for the cb coupling and of (9(A^) for the other two, A ~ 0.23 being 
the usual Cabibbo parameter. This conclusion also implies that current data allow deviations from 
unitarity in the quark-mixing matrix only at (9(A^) or higher. On the other hand, new physics 
affecting right-handed quarks can be of C(A^) for cd and cb transitions, and as large as 0{X) for 
cs transitions. That is to say that the right-handed cs matrix-element is the least constrained. 

Perhaps surprisingly, we note that the constraints displayed in Table II are comparable or tighter 
than existing constraints on anomalous VT-boson couplings to the top quark [3]. We can gain more 



19 
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FIG. 2: Parameter space of IniK^ and IniK^ subject to constraints from the contributions of magnetic- 
dipole operators to e, e', A'=\, and the neutron EDM for two representative sets of i?£^e',+,-,n- The various 
regions are described in the text. 



insight into these numbers by interpreting them in the context of left-right (LR) models with mixing 
of the and Wb. gauge bosons. In these models one predicts ^ ~^^w^ where is the VFl-W^r 
mixing angle. This angle is constrained by b ^ s'~f to be at the 10~^ level [40], and so in these 
models is only allowed at the 10~^ level. The additional freedom found in our study arises from the 
general decoupling between the top and charm anomalous couplings. 

For the right-handed anomalous charm couplings, the LR models result in the generic form 
nf^ ^ {g^i/ ghj^w^ch/^cD D = d,s,b. The first factor, {g^/ gi}^w^ allowed to be several 
times larger than [40], whereas the second factor depends on the right-handed mixing matrix, V^. 
Our bounds in Table II suggest within this context that constraints on right-handed mixing-matrix 
elements involving charm, V^, are not very tight at present, with V^^ being the least constrained one. 

Finally, our study also indicates which future measurements provide the most sensitive tests for 
new physics that can be parametrized with anomalous charm- couplings. For the CP- violating 
imaginary parts, the n EDM and the hyperon asymmetry are the most promising channels 
for probing right-handed couplings. To probe CP-violating left-handed couplings, more precise 
measurements of sin(2/5) and sin(2/3s) are desired. Constraints on the real parts of the right-handed 
couplings can be tightened with improved measurements of semileptonic B and D decays. 
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APPENDIX A: INDEPENDENT PARAMETERS IN QUARK MIXING MATRICES 



The parametrization we introduced in Eq. (1) reflects two effective matrices, K^^, for the charged 
currents involving left- and right-handed quarks, respectively. In the effective theory under consid- 
eration, these 3x3 matrices are complex and nonunitary; the usual unitarity of the CKM matrix is 
lost. The general 3x3 complex matrix has 18 parameters: 9 magnitudes and 9 phases. Nevertheless, 
not all the parameters in the matrix describing the left-handed charged current are independent. 
Of the 9 phases, 5 can be removed by redefinitions of the quark fields. Considering a scenario in 
which the corrections to the CKM picture in the SM are small, we find it convenient to choose these 
parameters as 



K, 



eff 



ud 



Kb e 



Vll 



vji 



V 



ts 



tb 



\ 



J 



(Al) 



More explicitly, four of the five quark phases have been used to remove the phases in Vud, Kis, 
Vtsi and Vtb- The remaining quark phase has to be used to remove one of the phases in the charm 
row, and for convenience we choose it to be 0^^, thus setting it to zero. Equivalently, only the two 
relative phases between the three (-K"^) elements are physical. These, plus the phases (p^^ and 
of \4b and Vtd which become independent, are the four physical phases in the general (nonunitary) 
3x3 complex matrix. In this study we only allow for three physical phases by requiring that the 
phases of Vub and Vtd be related, as in the usual CKM picture. 



There are also nine independent real parameters, the magnitudes of each of the elements of 



cflf- 



For our study we concentrate on the possibility of new physics in the charm sector, leading us to 
reduce the number of free parameters that we consider from 9 to 6 by simply setting the other 3 to 
zero. We choose 3 of the 6 to be the 3 angles that describe the usual CKM matrix and the rest the 
3 anomalous couplings that appear explicitly in K]^^. We are thus left with 



K. 



efT 



A2 
2 



-A(i + h^,|) 



1-^ 

2 



\A\^{1- p-lT]) 



A 



A\^{p - ir]) 



(A2) 



where, for simplicity of notation, we have included here only the usual terms up to order in 
the Wolfenstein parametrization of the CKM matrix; the complete, unitary CKM matrix should 
be understood as remaining in Eq. (A2) in the limit and (p^j 0. In this paper we 

have provided formulas for the relevant observables in terms of these nine parameters. This set 
of formulas would allow us, in principle, to repeat the global fits for this scenario and constrain 
the nine parameters. We content ourselves with a less ambitious analysis in which we assume, in 
accord with observation, that the CKM picture is dominant and deviations are small. We thus use 
the CKM parameters as extracted from the global fits of CKMfitter [28] as input. The subsequent 
comparison of this global fit with specific observables then reflects the extent to which deviations 
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from unitarity are allowed by current data. Conservatively, we carry out this comparison at the 
one-sigma level. 

None of the above considerations apply to the right-handed charged current: all the SM field 
phases are fixed by the removal of five phases in K^'f^ so that all 18 parameters in the corresponding 
are physical. In this study we limit ourselves to 6 of them, which we write as 



K. 



eff 



v 















-A K 



^1 e^^^ 















(A3) 



As pointed out in the conclusion, one needs to keep in mind with this parametrization that the new 
physics is not just hi-j but its product with V^j, as explicitly seen in Eq. (A3). 



APPENDIX B: AMPLITUDES IN UNITARY GAUGE 



The loop diagrams that are relevant to some of the processes we consider are displayed in 
Figs. 3 and 4. In evaluating the diagrams, we use dimensional regularization with a completely 
anticommuting 7^ and adopt the unitary gauge, which implies that they contain only fermions and 
W^-bosons. Since the theory with anomalous couplings is not renormalizable, some of the one-loop 
results are divergent. For these we adopt the prescription to drop the combination 2/(4 — D) —7^ + 
ln(47r) in the D-dimensional integral, 7^ being the Euler constant, and retain the accompanying 
logarithmic, ln(yU,/m^), part as well as other finite terms that depend on the mass of the quark 
in the loop. Moreover, we identify the renormalization scale /x with the scale A of the new physics 
parametrized by the anomalous couplings, fi = A. In the SM limit (k^'^ = O) , after the unitarity 
relation V*^,V^^ + V*^,V^^ + V^*^,V^^ = is imposed, our results are finite and reproduce those obtained 
in the literature in gauges [33]. 




FIG. 3: Diagrams contributing to amplitudes for d — > d'V*, with V being a neutral gauge boson. In all 
diagram figures, straight lines denote fermions and the loops contain W bosons besides fermions. 



K/x/x/x/v- 

(a) (b) 

FIG. 4: Box diagrams contributing to amplitudes for (a) dd' — > i'^i" or I'v and (a,b) dd' — > dd' . 
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We will present the details of our calculation elsewhere [41]. Here we provide the resulting 
effective Hamiltonians relevant to the loop-induced processes dealt with in this paper. 

The effective Hamiltonians for electromagnetic and chromomagnetic dipole operators involving 
down-type quarks d and d! ^ d are derived from the four diagrams in Fig. 3 with up-type quarks q 
in the loops and can be expressed as 

= E d'a^^{FlP^ + FlP^)dF^^ , (Bl) 

q=u,c,t 

-H^^d', = E d'<y^,{GlP^ + GlP^)tJG^\ (B2) 

* q=u,c,t 

where af"" = ^[-ff'^Y], 

n = {>im,, + \fm,)F^^{x^) + A^(l + «a<;m,Fo(g , 



G'n = (Ajm, + Ajm,,)GrK) + A,(l + z^^^j^m^ Go W , 



with 



— 2 

m 



F^^(x) = z z Inx 

° ^ ' 24(1 4(l-a;)4 

-20 + 31x-5x2 2x-3x2 
^"^"^^ " 12a 2(1 -xY ' 



(B3) 



(B4) 



A, = V;,,V^, , = A Jl + (1 + k];,) , Af = A^ ^f:,^^ , x^ = , (B5) 



(B6) 



gj^ — 2a; — + 3x^ Inx — 4 — x — 3x Inx 

^"^^ = 8(1 4(1 -x)4 ' ^0^"^^ = 4(1 -x)2 " 2(1 • ^ 

We note that the loop calculation for these operators yields finite results. We also note that 
different notations, D'q = —2Fq^, Eq = —2Gq^, F = Fq, and G = Gq, are sometimes used in 
the literature [21, 22]. 

The process dd' vv receives contributions from all the diagrams in Fig. 3 via rfc?' Z* ^ vv 
and from the box diagram in Fig. 4(a), the loop fermions in the latter being an up-type quark and 
a lepton. After summing over q = u,c,t and imposing the unitarity condition A^ + Ac + Af = 0, 
we find the effective Hamiltonian 

u -T/SM I njK /'"^o^ 
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where 
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(B9) 



UK 



a Gt 



dd'—^uu 



8 TT sin^ e 



w 



V ''V / 



+ 



a Gt. 



8 TT sin^ 9^ 



(4xg - 3) In 



A 



m 



w 
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^ r H ^ TTT In X 

8{x-l) 8(x-l)2 



= 2x-— — — \nx - AXJx) 
1 — X 



(Bll) 



In the expressions for H^d'-^up above, we have neglected the dependence on the mass of the loop 
lepton in the box diagram, but it is possible to generalize the formulas to include the dependence 
on tau-lepton mass [41]. 

The amplitude for dd' — > gets contributions from all the diagrams in Fig. 3 via dd' — > 

{Y,Z*) i+i- and the dia gram in Fig. 4(a), the loop fermions in the latter being an up-type 
quark and a neutrino. The resulting effective Hamiltonian is 
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'^dd'^e+l- 



dd'^e+e- ' 



(B12) 



where 
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(B13) 
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(B14) 
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with 



x{x — 4) 

8(x-l) 8(a;-l)2 



^ , , 18x^ - 163x3 + 259x2 - lOSx 24x^ - 6x^ - GSx^ + 50x . , 
^o(^) — ^"TTT 7T^ 1 z::r7 ttt; inx 



(B15) 



144(x- 1)3 
5x — 2x2 



72(x- 1)^ 



F(x) = -2x + ^ lnx + 4ro(x), Z(x) = 2x - 4x In x + 8Zo(x) . (B16) 

In '^^^'Zt+i- ^-bove, we have not displayed terms contributed by the magnetic {a^") parts of the 
dd' — *• 7*, amphtudes for convenience, but they can be found in Ref. [41] and do not contribute 
to the decay yU'^/i" which we consider. We note that the SM results in Eqs. (B9) and (B13) 

are in agreement with those found in the literature [21, 33]. 

From the two box-diagrams in Fig. 4 with quarks d and d' in the external legs and quarks q and 
q' in the loops, we derive the effective Hamiltonian 



n_j _ _ 'XySM _ I n/K _ 



"■dd'^dd' I ' ^-dd'^dd' ' 



(BIT) 



where 



'^dd'^dd' 



F '"'W / \2 



47^2 



(a2^o(^c) + A?^oN +2AA^o(^c,^0) d'l^P^dd'^^P^d , (BIS) 



l^dd'^dd' - i«^2 Z^IW'^WJ 



with 



+ 



167r2 
47r2 
47r2 



q,q 



A 



(6-2xJln--^ - B,{x^,Xg,) 

"I'W 



dYPLdd'laPLd 



q,q 



(6 - x^ - Xg,) In - ^2 (x,, X,,) 

"I'W 



dYPLdd'^Ad 



^ f ^d^'d d'P^d d P^d + ti^d''^^*d' d Pi^d d'P^d j , 



-3x1/ 



xy (4 — 8x + x^) In x xy (4 — 8y + in y 



4(l-x)(l-y) 4(|/-x)(l-x)2 4(x - ?/)(l - y)^ 



(B19) 



(B20) 



3, , 3X + 7/-X1/ 4x2-8x3 + x4)lnx 4y2_8y3 + y4^1n2/ 
2 {l-x){l-y) {y-x){l-xf [x - y)(l - yf 

n(rr„\ = ^(r\,A ^i^ + y-^y) I (4-x)Vlnx (4-t/)Vlny 
' 2^"^ + (1 - x)(l - y) + (y - x)(l - x)2 + (x - y)(l - y)^ ' 

xy-x-y-2 (4x - 2x2 + x^) Inx Uy - 2y'^ + y^) Iny 
I3s{x,y) = — r + — — + 



(B21) 



(l-x)(l-y) (?/-x)(l-x)2 (x-y)(l-y)2 
and «S'q(x) = limj^^a, Sq{x, y). The expression for 'H^^^^^, agrees with that in the literature [21, 33]. 
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APPENDIX C: MATRIX ELEMENTS AND PARAMETERS 



The matrix elements used in estimating the amphtudes for B J /ipK^rj^K are 



c L 

2 



(CI) 



where A*"^ = 3 is the number of colors and m and n are color indices. The decay constants above are 
= 416 MeV extracted from T{J/tp e+e") data [13] and /^^ = 420 MeV calculated in Ref. [42]. 
The form factors are functions of {p — q)^, with F(f^(m^^) = 0.45 and F^^ {m^') = 0.65 

from Ref. [43]. For meson masses, we use the values in Ref. [13]. 

The K, Bdi and Bs decay-constants are [13, 28] 

= 155.5 ±0.8, fB, = 191 ±15, /b. = 228 ± 17 , (C2) 

all in units of MeV. All of the following parameter values are obtained from Ref. [28] , the experimen- 
tal and theoretical errors given therein having been combined in quadrature. The QCD-correction 
factors in the K- and -B^^s-mixing amplitudes are [21, 28, 44]. 

r/^, = 1.46 ±0.22, r/^ = 0.551 ± 0.007 . (C3) 

The bag parameters used in the i^'-mixing amplitude are defined by 

'i 3(m^±mJ 

and similarly for Bb^ ^ in the B^ ^-mixing cases, where [28] 

Bj^ = 0.72 ±0.04, = 1.17 ±0.08, B^^ = 1.23 ± 0.06 . (C5) 

The charm and top masses used in the loop functions are [28, 45] 

Th^{m^) = 1.29 ±0.04, Th^{m^) = 165 ± 1 , (C6) 

both in units of GeV. For the CKM parameters, we adopt in the Wolfenstein parametrization the 
central values [28] 

A = 0.8116 , A = 0.22521 , p = 0.139 , f] = 0.341 . (C7) 

In our numerical estimates, we use only the central values of the parameters above, as their errors 
being no more than 20% are within the intrinsic uncertainty of our analysis. 
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